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We have performed large-scale Monte Carlo simulations on a model describing a (2+l)-dimensional 
array of dissipative Josephson junctions. We find three distinct stable quantum phases of the system. 
The most ordered state features long-range spatial ordering in the phase 9 of the superconducting 
order parameter, but temporal ordering only in spatial gradients AO, not in 9. Significantly, the most 
ordered state therefore does not have 3D XY ordering. Rather, it features 2D spin waves coexisting 
with temporally disordered phases 9. There is also an intermediate phase featuring quasi-long-range 
spatial order in 9 coexisting with a gas of instantons in A9. We briefly discuss possible experimental 
signatures of such a state, which may be viewed as a local metal and a global superconductor. The 
most disordered state has phase disorder in all spatio-temporal directions, and may be characterized 
as a gas of proliferated vortices coexisting with a gas of AS- instantons. The phase transitions between 
these phases are discussed. The transition from the most ordered state to the intermediate state 
is driven by proliferation of instantons in A9. The transition from the intermediate state to the 
most disordered state is driven by the proliferation of spatial point vortices in the background of 
a proliferated AS-instanton gas, and constitutes a Berezinskii-Kosterlitz-Thouless phase transition. 
The model also features a direct phase transition from the most ordered state to the most disordered 
state, and this transition is neither in the 2D XY nor in the 3D XY universality class. It comes 
about via a simultaneous proliferation of point vortices in two spatial dimensions and AS-instantons, 
with a complicated interplay between them. The results are compared to, and differ in a fundamental 
way from, the results that are found in dissipative quantum rotor systems. The difference originates 
with the difference in the values that the fundamental degrees of freedom can take in the latter 
systems compared to dissipative Josephson junction arrays. 

PACS numbers: 74.81.Fa,05.30.Rt, 74.40.Kb, 74.50. +r 



I. INTRODUCTION 

In general, dissipation suppresses quantum fluctua- 
tions and may support states of spontaneously broken 
symmetry. A remarkable consequence of this is the 
dissipation-driven quantum phase transition in a sin- 
gle resistively shunted Josephson junction in which the 
phase difference is localized in a minimum of the pe- 
riodic Josephson potential.^ In the parameter space of 
Josephson coupling and dissipation strength, this corre- 
sponds physically to a phase diagram with one metallic 
phase and one superconducting phase. While the behav- 
ior of a single dissipative Josephson junction is theoret- 
ically well understood, the picture is less complete for 
spatially extended systems. Other than the fully disor- 
dered phase and the fully ordered phase expected from 
the single-junction system, the phase diagram of arrays 
of dissipative Josephson junctions is co njectu red to host 
additional phases in both one^^and two 1 2 * 8 * 9 ! dimensions. 
These new, exotic phases can broadly be characterized by 
having various combinations of global and/or local phase 
fluctuations or order. 

Most of the analytical works on similar models have 
been based on mean-field analyses or perturbative renor- 
malization group arguments. Since these approaches are 
valid in a limited region of the parameter space, in par- 
ticular regions far away from phase transitions, a non- 
perturbative approach is of importance. Previous nu- 



merical work on models of dissipative Josephson junc- 
tions has mostly focused on lower-dimensional systems. 
The first Monte Carlo simulation of a single dissipative 
Josephson junction was presented in Ref. [TDJ where a 
fluctuation measure of the imaginary-time path of the 
phase difference was introduced to characterize the lo- 
calization transition. Improved and extended results for 
the same model were later reported in Ref. [TT] For one 
spatial dimension, Ref. H2I reported four physically dis- 
tinct phases for a dissipative Josephson junction chain. 
This simulation was performed on a dual model and not 
directly on the phase degrees of freedom. A model for a 
(2+l)-dimensional [(2+l)D] dissipative Josephson junc- 
tion arrays (JJA) has been treated numerically by Ref. 
ITU1 In essence, their results support the simplest sce- 
nario for a zero-temperature phase diagram) 1 ' 15 ' with one 
phase with and another without spatio-temporal order. 
This is also what was found in a large-scale Monte Carlo 
simulation on the dissipative (2+l)D XY quantum rotor 
modelP^ 

Finally, our investigations are also motivated by a 
rather different physical system which can be described 
by a closely related model. In Ref. [TT], a quantum 
XY model with bond dissipation in two spatial dimen- 
sions was used to describe quantum critical fluctuations 
in cuprate high-T c superconductors. The principal result 
of analytical work on this model is that the dissipation- 
driven quantum critical point is local, in the sense that 
the fluctuation spectrum is frequency dependent but mo- 
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mentum independent.^ Although the physical system we 
have in mind primarily is that of a Josephson junction 
array, we return to a discussion of the possibilities of local 
quantum criticality later in the paper 1^ 

The purpose of this paper is to numerically investi- 
gate the phase diagram of a specific model of a (2+1)- 
dimensional dissipative Josephson junction array. We 
pay special attention to the manifest anisotropy that ex- 
ists between the spatial and temporal dimensions. To 
be specific, the fluctuations of the quantum paths of the 
phase gradients will be explicitly characterized in terms 
of roughening transitions, allowing us to consider the 
(temporal) localization transition separately from the on- 
set of (spatial) phase coherence. In particular, we will 
identify a partially superconducting phase with spatial, 
but no temporal phase coherence. This corresponds to 
a dissipative JJA which may sustain a nonzero Joseph- 
son current, but where one nonetheless has voltage fluc- 
tuations over each junction. We investigate two phase 
transitions where the spatio-temporal aspects are well 
separated and can be characterized in terms of either a 
spatial vortex-antivortex unbinding, or proliferation of 
instanton-like defects. We also discuss a direct quan- 
tum phase transition from an ordered state to a disor- 
dered state involving simultaneous disordering in space 
and imaginary time. This corresponds to a quantum 
phase transition on a dissipative JJA where one tran- 
sitions from a state sustaining a Josephson current and 
allowing no voltage fluctuations to a normal state, but 
via an unusual quantum phase transition that is neither 
in the 2D XY nor 3D XY universality class. 



A. Model 

An array of Josephson junctions consists of supercon- 
ducting islands arranged in a regular network. Separat- 
ing the islands arc tunnel junctions in which Cooper pairs 
are able to tunnel from one superconducting grain to the 
neighboring grain. The fundamental degrees of freedom 
are the phases of the superconducting order parameters 
residing on the grains. A classical two-dimensional JJA 
is described by the 2D XY model 



H 



-K 



E ' 

<x,x'> 



i(6 x — X >), 



(1) 



where the summation goes over nearest neighboring sites 
on a square lattice. X is the phase of the complex order 
parameter of the superconducting grain at position x. 
Although the U(l) symmetry of the phase variables can- 
not be spontaneously broken in two dimensions at any 
nonzero temperature (implicit in the classical descrip- 
tion), the system nevertheless undergoes a Berezinskii- 
Kosterlitz-Thouless (BKT) transition in which it de- 
velops quasi-long-range order (QLRO) with power-law- 
decaying correlation functions in the low-temperature 
regime. The low-temperature phase corresponds to a 
dipole phase where the vortices and anti-vortices of the 



phase field are bound in pairs. At the transition the vor- 
tices proliferate and destroy the QLRO. For a given phase 
configuration, a single vortex is identified on a plaquette 
by a nontrivial line integral of the phase difference around 
the plaquette, taking the compactness of the phase field 
into account. 

The quantum generalized version of the model includes 
two additional terms describing quan tum fluctuatio ns in 
imaginary time r. The action read d 3 | 12 | 13 | 15 | 17 | 19 | 2 Sl 
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where we have defined the lattice gradient A# XX '. T = 
$x,r — Ox'.t- The first term describes the self-capacitance 
of a single island, the second term is the familiar Joseph- 
son interaction, coupling each superconducting island 
to the nearest neighbors by a periodic potential. The 
last term describes the ohmic dissipation as modeled 
by a bath of harmonic oscillators coupling to the bond 
variables J^S 

A subtle consequence of the presence of this ohmic 
shunt mechanism is that the phase variables become 
noncompactpZl as the dissipation term in Eq. Q breaks 
the 27r-periodicity of the Josephson potential. Thus, 
the phases are no longer defined with compact support 
9 6 [— 7r,7r), as they would be in the non-dissipative case 
or in a (2+l)D dissipative quantum rotor model. In- 
stead, we have 6 G (—00,00). The impact of this de- 
compactification on the problem is enormous. It reflects 
that a sudden increase along imaginary time in the phase 
difference, e.g., A9 X , X ', T — > A8 x ^' tT + 2n, would pro- 
duce a voltage imbalance over the barrier. A dissipative, 
measurable current would then flow through the shunt- 
ing resistors until the imbalance is relaxed. Hence, the 
variables cannot be defined modulo 2n, since A0 x x ' r 
and A0 X!X ' !r + 2-7T represent distinguishable states. The 
noncompactness of the variables implies that we may no 
longer identify vortices in the same manner as described 
above, as a line integral around a plaquette always yields 
zero for a noncompact phase field. In App. [X] we intro- 
duce a reformulation of the phase variables in terms of a 
compact part and an additional field describing the tun- 
neling between wells in the extended Josephson poten- 
tial. This enables us to identify vortices in the compact 
part of the phase. The phase transitions involving spatial 
ordering may therefore still be described by vortex pro- 
liferation even though the variables are of a noncompact 
nature. 

As a description of a dissipative JJA, there are a 
few simplifications built into the action We have 

only considered the effect of self-capacitance and ne- 
glected mutual capacitive coupling with neighboring 
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grains. Also, the dissipation term only accounts for one 
source of dissipation, namely the flow of normal elec- 
trons through the shunting resistors. Additional dissipa- 
tive effects like quasiparticle tunneling^ and Cooper pair 
relaxatiorP*22I have been neglected. 

In order to study the behavior of a two-dimensional 
array of Josephson junctions at zero temperature under 
the influence of ohmic dissipation, we perform large scale 
Monte Carlo simulations on a discretized version of Eq. 
d2k 



S =-y 1 ^ E^ X ' T+1 - O^tY 



N JV T 



(3) 



x r 



- K E X)cos(Ae x , x , )T ) 

<X,X'> T 



E E 

(x,x') t^t' 



7T V ( A ^x,x',r - A0 x , xV ,) 2 



,\N T J sin^ \t-t'\) 



Here, K T = 1/BpAr and the spatial coupling has been 
renamed KAt — > K. Our goal is to investigate the be- 
havior of the system in the K-a-spa.ce, K T therefore de- 
fines the energy scale and will be kept at suitable val- 
ues in the simulations. The variables are defined on the 
vertices of a three-dimensional cubic grid. The spatial 
linear extent of the grid is given by N, and the number 
of Trotter slices used to discretize the temporal direction 
is given by N T . Thus, At = f3/N T , and the size of the 
space-time lattice is N x N x N T . Periodic boundary 
conditions in imaginary time are implicit from the path 
integral construction, and are also applied in the spatial 
directions in the standard manner. The noncompactness 
of the variables also dictates the form of the kinetic term. 
Because 8 is an extended variable, its derivative must be 
expressed by discretized differentiation. We refer to the 
appendix of Ref. for details. 



B. Outline and overview of main results 

For outlining the roadmap to this paper, the phase dia- 
gram of the system is helpful. This is illustrated schemat- 
ically in Fig. 1. In all regions of the phase diagram, the 
phases 9 are disordered in the imaginary-time direction. 

In Sec. ^ we introduce the various observables used 
to identify the phases and phase transitions of the model 
defined in Eq. In Sec. Ill the details of the Monte 



Carlo simulations are presented in a concise form. 

In Sec. |IV[ we take a large value of the Josephson 
coupling K and investigate the behavior of the system 
as it crosses from the CSC phase to the FSC phase in 
Fig. [T| upon increasing a. There is a phase transition 

(2) 

at a critical dissipation strength, a c , above which the 
system is fully bond-ordered superconducting (FSC). For 

(2) 

a < a c the system features unbounded temporal fluc- 
tuations, while at the same time featuring spatial phase 
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FIG. 1: A schematic phase diagram of the system defined 
by Eq. [3] based on the Monte Carlo calculations presented 
below. Here, we have used the value K T = 0.002, corre- 
sponding to the parameters in Sec. |IV| NOR refers to the 
normal phase, where vortices are proliferated and the bonds 
AO are disordered in the r direction. CSC refers to the critical 
superconducting state, where the 9 variables feature power- 
law correlations in space, while A6 remains disordered. FSC 
refers to the fully bond-ordered superconducting state, which 
features an additional ordering compared to the CSC phase, 
namely A(9-ordering in imaginary-time direction. A hypothet- 
ical fourth, local phase has not been observed in our simula- 
tions, as indicated by the box in the lower right corner. See 
text in Sec. II Bl for more details. 



coherence. Due to algebraically decaying spatial correla- 
tions in this regime, we will refer to the phase as critical 
superconducting (CSC). In other words, the phase con- 
figurations of the system rotate more or less as a "rigid 
body" in time, thus at the same time giving rise to a fi- 
nite superfluid density (helicity modulus) as well as volt- 
age fluctuations across the junctions. A detection of the 
CSC phase thus requires simultaneous measurements of 
the superfluid density of the system, as well as ac mea- 
surements of voltages across junctions. 

In Sec. [V] we consider the transition between the NOR 
phase and the CSC phase, and this is found to be a purely 
spatial phase transition of the BKT type. 

In Sec. |VI[ we investigate the response of the system to 
increasing dissipation at low and intermediate Josephson 
couplings, as it crosses from the NOR phase to the FSC 
phase in Fig. [T] This is the most difficult case to ana- 
lyze, as the system transitions from a spatio-temporally 
disordered phase directly to the to the spatio-temporally 
ordered state FSC upon crossing the critical line etc ■ 

In Sec. |VII| the topological defects driving the various 
phase transitions as well as how such a model may exhibit 
local quantum criticality (LQC), are discussed. This may 
be briefly summarized as follows. 

On the line separating CSC from NOR, and on the 
line separating CSC from FSC, the spatial and temporal 
aspects of the phase transitions can be considered sepa- 
rately. The CSC-NOR transition is driven by point vor- 
tices and is in the 2D XY universality class. The FSC- 
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CSC transition is driven by instantons in A9 and may be 
characterized as a roughening transition in the space of 
AO. On the critical line there is a complicated in- 
terplay between temporal and spatial fluctuations. This 
critical line is neither in the 2D XY nor in the 3D XY 
universality class. 

An additional fourth phase could conceivably have 
been present in the phase diagram, featuring temporal 
order and unbound vortices. The most likely position in 
the phase diagram for such a hypothetical phase would 
be at weak Josephson coupling and strong dissipation 
strength. This is shown by the dotted lines within the 
box in lower right corner of Fig. ^ The local transition 
line would involve ordering of temporal fluctuations with- 
out onset of spatial phase coherence, and as such would 
describe a local quantum critical point. Our simulations, 
however, show no sign of such behavior in the parameter 
range we have considered. 

The limit a — is in principle ill-defined in this model 
since a finite dissipation is essential for the decompact- 
ification of the variables. This is indicated by drawing 
the a = axis as a red dotted line in Fig. [T] Thus, the 
value a = is also a singular endpoint of the horizontal 
(red) line in the phase diagram, and this is indicated by 
terminating this line in an arrow. 

In App. [X] we provide some more details and dis- 
cussion on the fundamental implications of the noncom- 
pactness of the phase field. In App. [Bj we take a closer 
look at the NOR phase and investigate the description 
of Refs. [23 E2 of such a normal phase as a so-called 
floating phase. 



II. OBSERVABLES 

In order to describe the various phases and transitions 
introduced in the previous section, several quantities will 
be calculated. To monitor the degree of (spatial) super- 
conducting order, we calculate the spatial hclicity mod- 
ulus, or phase stiffness. This quantity measures the in- 
crease in the free energy when applying an infinitesimal 
twist across the system, 9 X — >• 9 X — 5 ■ x. It probes the de- 
gree of phase coherence in the system and thus its ability 
to sustain a supercurrcnt. The only term in the action 
that contributes to the helicity modulus is the Joseph- 
son interaction term. Hence, the helicity modulus is 
given by 
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Here, the brackets indicate ensemble averaging. In the 
context of the classical 2D XY model, = defines 
the disordered state where vortices are proliferated. In 



the same manner, T x signals the finite rigidity of the 
quasi-ordered state. 

The same XY models used to describe superconduct- 
ing systems also describe magnetic systems of planar 
spins, and the superconducting phase 9 can formally be 
associated with the direction of the XY spins. Conven- 
tionally, the order of a superconducting system is there- 
fore often described by a magnetization order parameter 
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which probes the uniformity of the spin direction across 
the entire (2+1) -dimensional volume of the system. 

It should be noted that these two order parameters are 
periodic and consequently insensitive to tunneling events 
where the phase difference on a single junction jumps to 
a neighboring potential well, A9 — > A9 + 2ir. Conse- 
quently, T x and m do not probe the dissipation-induced 
localization per se. In order to quantify this, we calcu- 
late the mean square displacement (MSD) of the bond 
variable A9 along imaginary time, 
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Here, we have defined Ad = l/N T J2 T A0 T . The MSD 
is often used in the context of stochastically growing in- 
terfaces or diffusion processes, and it is natural to adopt 
some concepts from these areas for our problem. For in- 
stance, the degree to which the imaginary-time history of 
A9 may be regarded as "rough" can be quantified by the 
scaling characteristics of the MSD with the length N T of 
the "interface". Normally, one finds 
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if the imaginary-time history of A9 describes self-affine 
configurations. H = 1/2 corresponds to a Markovian 
random walk, and such linear scaling of the MSD is also 
referred to as normal diffusion. A deviation from linear 
growth of W\ 6 as a function of N T is the hallmark of 
anomalous diffusionW^In particular, H < 1/2 is referred 
to as subdiffusive behavior. A smooth interface is char- 
acterized by the MSD being independent of the system 
length. 

To describe the phases and phase transitions, we will 
also investigate correlations of the order parameter field 
considered in Eq. ([5| . We define the spatial and temporal 
correlation function by 
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where /i S {x,r}. The extra factor q in the exponent is 
introduced for later reference in App. [Bj but will be set 
to the conventional value q = 1 otherwise. In App. |Bj 
we will also consider bond correlations, defined here for 



convenience as 
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For completeness we also present the susceptibility of 
the action, 
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as an additional means of locating the expected 
dissipation-induced phase transitions. This is the quan- 
tum mechanical equivalent of the classical heat capacity 
and is expected to present a nonanalyticity at a critical 
point. 



III. DETAILS OF THE MONTE CARLO 
CALCULATIONS 

Considerable progress has been made in construct- 
ing new, effective, non-local algorithms fo r long-r ange- 
interacting systems with extended variables! 11 * 27 * 28 ^ How- 
ever, these algorithms are presently restricted to (0+l)D 
systems, and do not seem to generalize easily to N > 
1^3. In the Monte Carlo simulations, we have there- 
fore com bined local updates with a parallel tempering 
algorith m 129 ! 30 ! in which several systems are simulated si- 
multaneously at different coupling strengths. 

A Monte Carlo sweep corresponds to proposing a lo- 
cal update by the Metropolis-Hastings algorithm sequen- 
tially for every grid point in the system. The proposed 
new phases are generated by first randomly choosing to 
increase or decrease the value, then propagating the value 
by a random increment of size 27rn/32, where n £ {1, 32}. 
In other words, the continuous symmetry of the variables 
is emulated by 32 discrete states per 2ir interval. We have 
confirmed that adding additional states will not change 
the results. 

After a fixed number of Monte Carlo sweeps a paral- 
lel tempering move is made. In this move, a swap of 
configurations between two neighboring coupling values 
is proposed, and the swap is accepted with probability 
E PT given by 



1 if A < 0, 
- A if A > 0. 
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Here, A = k'(S[X;k']-S[X';k'])-k(S[X;k]-S[X';k]), 
where k is the coupling value varied, representing in our 
case a or K, and X represents the phase configuration. 
S indicates the part of the action conjugate to the cou- 
pling parameter k. Both the Metropolis updates and the 
parallel tempering swaps are ergodic and respect detailed 
balance. 

All Monte Carlo simulations were initiated with a ran- 
dom configuration. Depending on system size, various 
numbers of sweeps were performed for each coupling 
value. Error bars are provided for all observables except 
correlation functions, but are usually smaller than the 
data points. Measurements on which we perform scaling 
have, broadly speaking, a relative error well below 1%. 
The Mersenne- Twister 3 -^ random number generator was 



used in all simulations and the random number genera- 
tor on each CPU was independently seeded. It was con- 
firmed that other random number generators yielded con- 
sistent results. In some simulations we also made use of 
the Ferrenberg-Swendsen rewcighting technique, 32 which 
enables us to continuously vary the coupling parameter 
after the simulations have been performed. 

In order to identify sharply defined nonanalyticitics 
and observe converged scaling of W^ g at the dissipation- 
induced phase transitions, relatively large values of N T 
are needed. This limits the range of spatial sizes acces- 
sible in simulations with a single-site update algorithm. 
In the sections where we focus on the temporal scaling, 
we have fixed the spatial size at N = 20 and varied the 
temporal size in the range N T = 50 to N T — 350. In Sec. 
[V]we find that in the CSC phase the temporal size of the 
system is irrelevant in determining the spatial properties 
of the system. Consequently, the temporal size is fixed 
at N T = 20 and the spatial size is varied in the range 
N = 10 to N = 100. To investigate the spatial corre- 
lations in 9 across the NOR-FSC phase transition, we 
have also performed simulations on a N T = 30 system 
with N = 50 and N = 100. 



IV. THE CSC-FSC TRANSITION a< 2) 

In this section, we consider the behavior of the system 
under the influence of strong Josephson coupling K, i.e. 
for a value of K where the corresponding classical sys- 
tem would be topologically ordered even in the absence 
of dissipation (above the horisontal line in the phase dia- 
gram of Fig. [T]). The coupling parameter will be fixed at 
K = 1.5 in this section, while the dissipation strength a 
is varied. We will use a quantum coupling K T = 0.002. 
The main focus is on scaling of observables describing 
temporal fluctuations. Hence, the spatial system size is 
fixed at N = 20. 

We start by presenting typical configurations of the 
bond variable, A9, as a function of r. At strong cou- 
pling, the bond variables are located predominantly in 
the vicinity of the potential minima located at 2irn, 
where n is an integer. Due to the noncompact nature 
of the variables, A9 are free to tunnel between neigh- 
boring minima at weak dissipation. When considering 
the single-junction problem, this sudden tunneling of 
the bond variable from one Trotter slice to the next, 
A# T+ i — A9 T w 27rni, is often referred to as instanton 
or anti-instanton configurations, depending on the sign 
of the integer valued "instanton charge", n\. Note that 
the noncompactness allows for tunneling of A9 also be- 
tween minima of the potential located further away than 
nearest neighbor. This corresponds to instanton charges 
with values larger than unity. The tunneling behavior is 
easily identified in the topmost curve in Fig. [2] where 
frequent instantons and anti-instantons are apparent. In 
this temporally disordered state, the quantum paths of 
A9 appear to be well described in terms of a gas of pro- 



6 



8 




-1 I ' ' ' ' ' ' ' ' ' 1 

20 40 60 80 100 120 140 160 180 200 

T 



FIG. 2: The bond variable A9 as a function of imaginary time 
r for two different values of dissipation strength, a = 0.0102 
and a = 0.0281, in the strong Josephson coupling regime. In 
the topmost curve, the bond variable clearly spends most of 
the time in the vicinity of the potential minima, although tun- 
neling events between minima are frequent. The lowermost 
curve corresponds to the fully bond-ordered superconduct- 
ing state where AG is localized and W\ e does not scale with 
N T . Note that in the CSC phase (topmost curve) the quan- 
tum paths of A8 are well described in terms of instantons 
where the fluctuations in imaginary time are mostly given by 
integer multiples of 2n, in contrast to the situation for the 
corresponding quantum paths of A8 in the NOR phase, see 
Fig. © 

liferated instantons. Beyond a threshold value of a, we 
observe a localization of A 8 in one of the minima of the 
Josephson potential. The imaginary-time history of a 
bond variable corresponding to this phase forms an es- 
sentially smooth surface and is given in the lower curve 
in Fig. [2j However, even though the phase gradients 
are localized, closely bound pairs of instantons and anti- 
instantons may still be present. 

In Fig. [3j we show the mean square displacement as 
a function of dissipation strength. The temporal bond 
fluctuations are clearly suppressed for increasing values 
of a. The different curves represent different values of 
N T . Two regions of different scaling behavior of W\ e 
as a function of N T can immediately be discerned. For 
weak dissipation, W\ B increases with N T , while W\ e is 
independent of the temporal size at strong dissipation. 
Separating the two regions is a precipitous drop in W\ B 
at a value of a that we will identify as the localization 
transition point ai 2 K 

Further information on the delocalized phase (CSC) 
can be found from investigating the dependence of W\ g 
on the temporal system size N T . Here, the MSD scales 
with 7V T according to 

Wl e = a{a)\nN T , (12) 

where a(a) is a continuously varying proportionality con- 
stant. In Fig. [4] we have plotted W^ 9 as a function of 
lniV r for various dissipation strengths. All but the low- 
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FIG. 3: (Color online) W^g, Eq. Q, as a function of dis- 
sipation strength for a system with K T = 0.002, K = 1.5, 
N = 20 and various values of N T . Note the kink in the curves 

at ot = ct c and the saturation of W Ag at a finite value for 

f 2) 

a > a c . Error bars are smaller than the data points. Inset: 
A blow-up of the region around ai . 



est curve represent dissipation strengths well below a c . 
A clear logarithmic dependence is seen for all values of 
dissipation strength in the CSC phase. The lowest curve 
with zero slope corresponds to a > a c , where temporal 
fluctuations are effectively quenched and does not 

scale with N T . In this way the increase of temporal fluc- 
tuations in A0 T as a is lowered may also be interpreted 
as a roughening transition at which the profile described 
by AO changes from smooth to rough. However, it should 
be noted that the logarithmic scaling presented in Fig. 
[3] does not conform to the scaling ansatz ([T]) for a self- 
affinc interface. Instead, AO is anomalously diffusive in 
the sense that H = 0. This is sometimes referred to as 
superslow diffusiorJ^In comparison, Ref. 1 10 1 found that 
for the corresponding normal phase of a single resistively 
shunted Josephson junction, the MSD follows the power 
law ([7]) with the exponent decreasing continuously with 
dissipation strength (from H = 1/2 for a > to H sa 
for a = a c ). 

As we show in Fig. [5] W/\e as well as the action suscep- 
tibility xs an d the helicity modulus all feature nonan- 
alytic behavior at the critical value ai 2 \ Fig. [5] therefore 
supports the notion that the transition at etc is indeed 
a genuine dissipation-induced quantum phase transition. 
Since we have shown that the system for a > a c has 
both spatial phase coherence and temporal localization 
of AO, we can identify this region as a fully bond-ordered 
superconducting (FSC) phase. However, T x > even 
(2) 

for a < a c , indicating that also the weak-dissipation 
CSC phase features spatial phase coherence. The kink 
in the helicity modulus shown in the inset in Fig. [5] may 
be attributed to the (slightly) reduced spatial rigidity 
as the bond variables delocalize in imaginary time when 
leaving the FSC phase. An important conclusion to be 
drawn from this, is that in the regime of strong Joseph- 
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FIG. 4: W^g as a function of In N T for various values of the 
dissipation strength a ranging from the weak-dissipation limit 
to the ordered state at the top and bottom, respectively. The 
dotted lines indicate the logarithmic growth of W\ e . Error 
bars are much smaller than the data points. 
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FIG. 6: (Color online) The spatial correlation functions 
G e (x- q = 1), Eq. @, calculated for a = 0.005, K T = 0.002, 
N T = 20, and for two values of the Josephson coupling and 
two values of the spatial extent TV. Both coupling values corre- 
spond to the CSC phase. The dotted lines show the power-law 
fit of the correlation functions. 
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FIG. 5: (Color online) The spatial helic ity modulus T x , 
Eq. Q, action susceptibility \s, Eq. (10 1, and mean fluc- 
tuation Wao, Eq. as a function of dissipation strength 
a for a system in the strong Josephson-coupling regime with 
K T = 0.002, K = 1.5, TV = 20 and TV T = 250. Inset: A blow- 
up of the helicity modulus around the a c transition. For 
a > a c , the dissipation renormalizes the spatial coupling 
strength so that a kink in T x is visible at the localization 
transition. However, this renormalization is miniscule. Pro- 
liferation of instantons across the line a[ 2 ' in Fig. [ljdoes not 
trigger a proliferation of vortices. 



son coupling, proliferation of instantons does not trigger 







(2) 

a proliferation of vortices at a c in Fig 

A possible physical interpretation of the behavior at 
strong Josephson coupling and weak dissipation is a 
phase where there are fluctuations of voltage (and thus 
also of normal currents through the shunts) even though 
a finite superfluid density allows the system as a whole 
to sustain an unimpeded supercurrent. For reasons that 
will be apparent in the next section, we have chosen to 



refer to this state as a critical superconducting (CSC) 
phase. Similar conclusions have been made earlier for 
(1+1)D systems, e.g., in Refs. IHEl andlT^l where the au- 
thors claimed to have found an additional superconduct- 
ing state characterized by spatial coherence but large lo- 
cal fluctuations. An experimental signature of the FSC- 
CSC phase transition would be to measure an abrupt in- 
crease in voltage fluctuations across each junction while 
the system maintains a Josephson current across the sys- 
tem as the dissipation strength is reduced. The phase 
CSC therefore represents a locally metallic (on each junc- 
tion) and globally superconducting (throughout the sys- 
tem) state. 



V. THE NOR-CSC TRANSITION 

We next consider the phase transition separating the 
CSC phase from the fully disordered state NOR. First, we 
note that the region at weak dissipation and low Joseph- 
son coupling in the phase diagram of Fig. [T] is spatially 
phase incoherent, T x — 0. This is therefore identified as 
the normal, metallic phase (NOR) of the dissipative JJA. 
To verify that the CSC state identified in the previous 
section by its finite spatial coherence is a distinct phase, 
we next show that it is separated from the NOR phase 
by a genuine phase transition and not just a crossover 
caused by the limited spatial extent of the systems. 

In Fig. [6j we show algebraically decaying correlation 
functions in the spatial direction in the CSC phase, indi- 
cating QLRO within each Trotter slice. In combination 
with the observation of vanishing order in the temporal 
direction (as measured by W^ g ), this motivates an inter- 
pretation of the CSC phase as a dimcnsionally reduced 
critical phase in which the Trotter slices are decoupled 
from each other. We verified that varying N T had no 
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impact on the results for any of the observables probing 
spatial behavior. Thus, the extent of the systems is fixed 
at N T = 20 in the following. 

We anticipate the phase transition separating the NOR 
phase from the CSC phase to be in the BKT universal- 
ity class. At the transition point, the helicity modulus 
is expected to scale according to the finite-size scaling 
function^ 



T X (N) = T x (oo) 1 



1 



1 



2 In N + C 



(13) 



where T x (oo) is the value of the helicity modulus as 
N — > oo and C is an undetermined constant. The critical 
value K c may be extracted by varying K until an opti- 
mal fit is achieved. In addition, at a BKT transition, the 
value of T^(oo) obtained at optimal fit should satisfy the 
universal relation Y x (oo)K c = 2/it. 

By treating both parameters as variables in the fitting 
procedure, no a priori assumption on the value of the 
jump is made. This value may consequently be used as 
an additional check on the validity of the conjecture of 
identifying the transition as a BKT transition. 

In Fig. [7] we present for various spatial system 



7(a) 



sizes and the corresponding fit with Eq. (13 1. Fi, 
shows results for the dissipationless limit, a = 0, 35 while 
Fig. |7(b)| gives the corresponding results for a — 0.005. 
At both dissipation strengths we observe optimal fit at 
K w 1.12. The insets presented in both figures show 
r T x {po)K, which should be compared to the broken line 
indicating the expected 2/ it universal jump of a BKT 
transition. These results demonstrate that the NOR 
CSC transition is a BKT transition. The temporal in- 
teraction terms are evidently completely incapable of es- 
tablishing temporal order at this transition. In partic- 
ular, when comparing Fig. 7(a) and Fig. |7(b)| corre- 
sponding to no dissipation and weak dissipation, respec- 
tively, no significant difference is visible. Even though 
the dissipation term has a major impact on the temporal 
fluctuations,^ the spatial helicity modulus appears com- 
pletely unaffected by the presence of dissipation in the 
CSC phase. 



The classification of the NOR-CSC transition is im- 
portant in two respects. Firstly, the finite-size analysis 
shows that the existence of a finite helicity modulus in 
the CSC phase is not a mere finite-size effect. Secondly, 
the analysis places the transition in the BKT universal- 
ity class. This would not have been possible if there were 
a divergent correlation length in the temporal direction. 
Such an effect would have been likely to show up as a 
break-down of the scaling procedure. In this way the 
analysis gives an indirect verification that the transition 
is of a purely spatial nature, and that the CSC phase is 
temporally disordered and spatially quasi ordered. 
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(a)The spatial helicity modulus T x as a function of spatial 
system sizes N for a = 0.0, K T = 0.002 and various 
Josephson coupling values. 
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(b) Spatial helicity modulus T x as a function of spatial system 
sizes N for a = 0.005, K T = 0.002 and various Josephson 
coupling values. 

FIG. 7: Comparison of calculated values of t he spatial helicity 
modulus T x with the scaling function ( 13 1 for two different 



dissipation strengths. For both values of a, a good fit is ob- 
served at K w 1.12. Inset: The universal jump of the helicity 
modulus is expected to be 2/tt for a BKT transition. This 
value is indicated by a broken line in the insets. The univer- 
sal jump as calculated from the fitting procedure is shown to 
be in good correspondence with the BKT scenario. 



VI. THE NOR-FSC TRANSITION 

The transition line ai 1 ^ is the only transition line in the 
phase diagram that exhibits a simultaneous temporal and 
spatial order-disorder transition. Hence, it involves an 
interplay between instantons (or instanton-like objects) 
and vortices, but in a complicated way that is not easy 
to disentangle. 

In this section, the Josephson coupling strength will 
be fixed at an intermediate value of K = 0.4, for which 
a classical counterpart of our model would be well inside 
the disordered phase (T^ = 0). The quantum coupling 
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FIG. 8: The bond variable A9 as a function of imaginary 
time t for two different values of the dissipation strength, a = 
0.011 and ot = 0.021, in the weak and intermediate Josephson 
coupling regime. These values correspond to the normal phase 
and the ordered phase, respectively. The quantum paths of 
A9 in the normal phase (relatively low values of K) exhibit 
fairly slow variations in time, and are not necessarily well 
described in terms of instantons. Note the contrast to the 
quantum paths in the topmost curve above and the topmost 
curve in Fig. [2] 



is set to K T =0.1. Note that this differs from the value 
of K T used to compute the phase diagram in Fig. [T] 

Fig. [8] shows typical configurations of the bond vari- 
able AO as a function of r for two dissipation strengths 
corresponding to regimes where the model behaves quan- 
titatively different. The topmost curve corresponds to 
weak dissipation, with anomalous diffusive behavior of 
the value of A9. The lowest curve represents the regime 
of strong dissipation, where the imaginary time-history of 
AO is qualitatively less rough and where we can therefore 
show below that the bond variable is localized. 

We see from Fig. [9jthat the amplitude of the temporal 
bond fluctuations are rapidly decreasing with increasing 
a. At a critical value of the dissipation strength a — a' 1 ' , 
the MSD features a steep drop marking the localization 
transition where the tunneling of A9 is suppressed suffi- 
ciently to give the bond variables a well-defined value in 
imaginary time W* We are once again able to distinguish 
between two separate states based on the scaling prop- 
erties of the MSD. In Fig. [TUJ we present a plot of the 
MSD as a function of In N T for several values of a. The 
lowermost curve in the figure again represents the FSC 
phase, a > a£ , where the MSD is independent of N T . 
All other curves represent dissipation strengths below the 
localization transition, and for these a clear logarithmic 
scaling is observed. In this way, there are distinct delo- 
calized and localized regimes for the bond variable also at 
weak Josephson coupling, and the temporal fluctuations 
in each of them behave in exactly the same way as for 
strong Josephson coupling.^ 

To confirm that the temporal transition at a = etc 



FIG. 9: (Color online) W\ g , Eq. |6}, as a function of dis- 
sipation strength a for a system with K T = 0.1, K = 0.4, 
N = 20 and various values of N T . Note the kink in the curves 
at a = cii and the saturation of W^ e (N T ) at a finite value 
for q > a^ 1 '. Error bars are smaller than the data points. 
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FIG. 10: Wl e , Eq. (6}, as a function of lniV T for various 
values of a ranging from the weak-dissipation limit to the 
ordered state at the top and bottom, respectively. The log- 
arithmic behavior found at large N T is indicated by dotted 
lines. Error bars are much smaller than the data points. 



Fig. [TTJthc hclicity modulus T x . Note the abrupt man- 
ner in which the phase stiffness attains a finite value at 
a = aj; 1 '. Even though the spatial extent of the system 
is relatively small, there is no weak-dissipation tail which 
would have been visible for too small system sizes. In 
the same figure we also show the root mean square dis- 
placement, Wj^g, and the action susceptibility, xs- It is 
clear that all observables feature a nonanalyticity at the 
same point. We can therefore conclude that the transi- 
tion NOR-FSC is a quantum phase transition involving 
simultaneous onset of spatial and temporal order. 
In Fig. 



11 we note that the non-analyticity in on 



the line is brought out very sharply at the system 
also marks the onset of spatial ordering, we show in sizes we consider in this case, namely 20 x 20 x 250. 
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FIG. 11: (Color online) Spatial helicity modulus T x , Eq. Q, 
action susceptibility xs, Eq. ( 10 1, and mean fluctuation Was, 
Eq. |6} , as function of dissipation strengths for a system with 
K T = 0.1, K = 0.4, N = 20, and iV T = 250. Note that T x 
vanishes continuously at ai (no jump). 



Assuming hyperscaling and two diverging length scales £ 
(spatial) and £ r (temporal), we may write 



T, 



2-d- 



(14) 



Here we have introduced the dynamical critical exponent 
z defined by £ r ~ £ z . The sharpness can thus be ex- 
plained by the large system size and diverging length 
scale in the r direction. Very little finite-size effects may 
then be expected due to the limited spatial extent of the 
system, since d = 2 and the spatial correlation length 
drops out of the scaling. 

Ordinarily, it would have been natural to attempt a 
scaling analysis of this phase transition based on the 
Binder ratio 



Q 



(H 4 ) 
(H 2 ) 2 ' 



(15) 



in order to extract the dynamical critical exponent of the 
system, z. Here m is the magnetization order parameter 
of the superconducting phases defined in Eq. §5§. An 
ordinary quantum critical point is characterized by di- 
verging lengths in space and time, £ and £ T , respectively. 
The Binder ratio is then expected to scale according to 



Q = Q 



N N T 

7' 7" 



(16) 



The correlation lengths entering here are correlation 
lengths of the phase-correlation function, measuring 
correlations in the spatial and r directions. Thus, it 
should be possible to collapse the Binder ratio curves, at 
criticality, as a function of N T /N Z for the correct value 
of z. 

We have attempted such an analysis in this case, and 
failed. In our computations, we have been able to identify 
a diverging length scale £ based on the above scaling 
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FIG. 12: (Color online) Double-logarithmic plot of the spatial 
correlation functions Gs(x;q = 1) at values of a correspond- 
ing to above, close to, and below the NOR-FSC phase transi- 
tion. The relevant coupling values are K = 0.4 and K T = 0.1, 
and the system size is given by N = 100 and N T =30. At 
a ~ otc the correlation function appear to be linear at long 
distances indicating scale invariant spatial fluctuations. The 
dotted line indicates this linear behavior. Thus, £ — > oo, or 
at the very least £ > N/2. 



approach, but not a diverging length scale £ r . The reason 
is that in our model, the coupling in spatial directions is 
effective in ordering the phases 9, while the coupling in 
the t direction is only effective in ordering bond variables 
AO, while the variables never order in the r direction. 
One may therefore define a diverging length £ entering 
Eq. |16[ but not a diverging length £ T . A diverging length 
scale in the r direction may very well exist for the bond 
variables AO, but not for the phase variables 0. 

The onset of long-range order in the variables in the 
spatial directions may be described by the spatial cor- 



relation function Gg(x;q — 1), Eq. (|8|. In Fig. 12 
we present spatial correlations corresponding to dissipa- 
tion strengths slightly below the NOR-FSC transition 

{a < ai 1 ^), close to the transition (a w ), and slightly 
above the transition (a > a£ ). The spatial correlation 
length appears to behave as expected for a second order 
phase transition into a phase with long-range (spatial) or- 
der, implying that the NOR-FSC transition is associated 
with a diverging length scale in the spatial directions. 



VII. DISCUSSION 

Since the work of HertzJ^ quantum critical points are 
commonly characterized by their dynamical critical ex- 
ponent z. Underlying Hertz' scaling theory is Landau's 
notion that all relevant fluctuations of a system may be 
ascribed to fluctuations of an order parameter PSl This is 
evident when considering that the exponent z is defined 
from a divergence of a length scale of the order param- 
eter correlation function. Such a characterization may 
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therefore be insufficient when the critical point cannot be 
well described by one single order parameter, a problem 
which has been pointe d out in different cases in recent 
theoretical worksP^^l 

The model studied in this paper may be related to 
a problem of this kind in the sense that we are unable 
to find one single order parameter adequately describing 
the spatial, temporal, and spatio-temporal phase tran- 
sitions separating the NOR, FSC, and CSC phases in 
Fig. [T] To substantiate this, we show in App. [A] that 
the noncompact 9 variables may instead be formulated 
by a combination of a compact phase field 9 G [— 7r,7r) 
and an additional integer valued field k containing infor- 
mation on what 2ir interval the original variable belongs 
to. Using the reformulation of the 9 variables described 
in Appendix ^ it is clear that the magnetization order 
parameter m only probes the order of the compactified 
part of the phase, 9, but is completely oblivious to the 
state of the integer- valued field k. Since the state of this 
field describes whether or not the phase differences A9 
are localized, m is fundamentally incapable of describing 
the localization transition concurring with the onset of 
coherence of 9. As a result, we are unable to define a 
dynamical critical exponent z. 

The phase transition from CSC to FSC is primarily 
temporal in the sense that it only involves condensation 
of instantons from a state where the spatial topological 
defects are already tightly bound. However, this localiza- 
tion of A9 also contributes to spatial ordering by coupling 
the Trotter slices along imaginary time, thereby reduc- 
ing spatial fluctuations sufficiently to render the system 
behavior 3D. Accordingly, CSC-FSC is also of a mixed 
character, as the transition separates a phase with spa- 
tial QLRO (CSC) from a phase where spatial long-range 
order is established (FSC). 

The phase transition from NOR to CSC is of a purely 
spatial nature. As one increases the Josephson coupling 
for weak dissipation, this transition involves only the 
binding of the (spatial) vortex degrees of freedom while 
the (temporal) instantons remain proliferated. This con- 
clusion is supported by the signatures of a BKT-type 
transition found in Sec. [V] In this way, the system be- 
haves as a stack of decoupled two-dimensional layers in 
the CSC phase, each exhibiting critical fluctuations in 
the 9 field. 

The phase transition from FSC to NOR is much more 
complicated than the ones from FSC to CSC and from 
CSC to NOR, and appears to be of a type not pre- 
viously considered in connection with superconductor- 
metal phase transitions. Since one cannot characterize 
the anisotropy of the phase transitions quantitatively in 
terms of an exponent z, we resort to more qualitative 
considerations of the spatial and temporal degrees of free- 
dom. In the case of intermediate coupling, one has a 
concomitant binding of vortices and localization of A9 
upon entering the FSC phase from the NOR phase. This 
corresponds to the ordering of the degrees of freedom 
relevant to space (9) and time (k), respectively. Due to 



this simultaneity, we characterize the NOR-FSC phase 
transition as a mixed spatio-temporal phase transition. 
It is an interplay between two distinct types of topolog- 
ical defects (point-vortices and temporal fluctuations in 
A9) that determines the character of the phase transi- 
tion. This phase transition is therefore neither of the 
BKT type, nor in the 3D XY universality class. The for- 
mer is characterized by proliferation of point-like vortices 
in two dimensions, while the latter is characterized by 
the proliferation of (2+l)-dimensional vortex loops^JttS. 
Dissipation, and the associated decompactification of the 
9 variables, leads to a disordering of the 9 variables in the 
imaginary-time direction in all regions of the phase dia- 
gram. Decompactification essentially chops up the vor- 
tex loops into spatial point vortices and instanton-likc 
objects in A9, thereby destroying the Lorentz-invariant 
physics of vortex-loop proliferation at the quantum phase 
transition. 

In order to exhaust all combinations of spatial and 
temporal order/disorder, one could also imagine a fourth 
phase exhibiting temporal order without accompany- 
ing spatial phase coherence. This would correspond to 
W^ g = const, and T x = 0, i.e., a phase with localized 
bond variables and proliferated vortices. The most prob- 
able location of such a phase would be at weak spatial 
coupling and large dissipation strength, corresponding to 
the lower right corner of Fig. [l] This scenario opens the 
possibility of a purely temporal ordering coinciding with 
exponentially decaying spatial correlations upon entering 
this hypothetical phase from the NOR phase. Due to this 
locality, such a transition could be a possible realization 
of a local quantum critical point ("z = oo") in a spatially 
extended system. In order to emphasize that the exis- 
tence of this local phase is only a possibility that we have 
not actually found in our computations, we have drawn 
a box of solid lines around the specific region in Fig. [T] 
and indicated possible realizations of the phase transi- 
tions by dotted lines. Although the existence of su ch a 
phase has been conjectured by analytical wortP^l anc j 
there is numerical work supporting this viewj^ W e find 
no signatures pointing to the existence of such a local 
phase in any of the parameter sets considered. Rather, 
our results strongly indicate that a spatial coupling is 
always rendered relevant by a large enough dissipation 
parameter In this way, the localization of A9 will al- 
ways induce an onset of spatial phase coherence. This is 
equivalent to saying that instanton-like excitations will 
always proliferate prior to, or simultaneously with, the 
unbinding of vortices as the strength of dissipation, a, is 
reduced. Local quantum criticality (in the sense of hav- 
ing temporal critical fluctuations coinciding with spatial 
disorder) would follow from vortices proliferating prior to 
instantons as the disordered state (NOR) is approached 
from the fully bond-ordered superdonducting state (FSC) 
by reducing a. 

Finally, we compare the phase diagram found in this 
paper with the phase diagram calculated for a mode l sim- 
ilar to Eq. ([3]) using compact variables.^ Fig. 13 shows 
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FIG. 13: Comparison of the phase diagrams of the noncom- 
pact (topmost) and compact (lowermost) models. Topmost 
diagram: The phase diagram found in this work. All phases 
feature disordered 9 variables in the imaginary-time direction. 
A notable feature is the phase CSC where bound vortex anti- 
vortex pairs coexist with disordered bond variables A8 in the 
t direction. This is a consequence of the 6 variables being 
defined with noncompact support. This is only true for finite 
a as the quadratic form of the dissipation term is the origin 
of the decompactification. Consequently, the physics found at 
finite a cannot be analytically connected to the limit a = 0. 
The description at a = would require compact phases and 
thereby a loss of the instanton degrees of freedom. For a = 
there is a phase transition at a critical value of K, but this 
phase transition is in the 3D XY universality class, as in 
the lowermost diagram, a = is therefore a singular end- 
point of the horizontal line in the topmost diagram. Lower- 
most diagram: The phase diagram found for a (2+l)D bond- 
dissipative quantum rotor model with compact variables. In 
this case the diagram features only a single transition line 
where the the system undergoes a spatio-temporally isotropic 
(z = 1) phase transition in the 3D XY universality class. See 
Ref. [If)] for details. 



two schematic phase diagrams, and the following discus- 
sion pertains to their topology. The topmost diagram 
summarizes the results found in this paper, while the 
lowermost diagram is the phase diagram for the (2+ 1)D 
dissipative quantum rotor model. In the latter case, the 
diagram features one single phase transition line separat- 
ing a completely ordered state from a disordered state. 
The phase transition separating them is driven by a pro- 
liferation of vortex loops. This transition line is isotropic 
in space-time (z = 1) meaning that the entire line is in 
the 3DXY universality class. From the lowermost phase 
diagram it is clear that if we start in the limit of no dis- 
sipation, a = 0, and increase a for K > K^dxy, the 
dissipation term only contributes to further dampening 
the innocuous three-dimensional spin-wave excitations. 
This can only increase the superfiuid density in the or- 
dered phase. However, in the noncompact model the 
regime K > A'bkt, and weak dissipation, represents a 
phase involving both two-dimensional spatial spin waves 
and a proliferated instanton gas. Increasing a from this 
regime may therefore drive a phase transition because 
the dissipation term is effective in binding these tem- 
poral defects. Therefore, the feature of the phase dia- 
gram of the noncompact model that really sets it apart 
from the phase diagram of the dissipative 3D XY model 
(i.e., the compact case), is the existence of a phase at 
weak dissipation involving spatial ordering concomitant 
with temporal disorder. The resulting phase CSC has no 
counterpart in the dissipative 3D XY model, since in the 
latter model the phases 9 are compact. Compact phases 
9 promote vortex loops as the critical fluctuations, while 
noncompact phases 9 promote vortices and instantons as 
relevant fluctuations driving the phase transitions. 

The phase CSC corresponds to a resistively shunted 
Joscphson junction array which may sustain a finite 
Josephson current through the array, but nonetheless fea- 
tures finite voltage fluctuations across each junction of 
the junction array. This may be viewed locally (at a sin- 
gle junction) as a metallic state, but globally (throughout 
the system) as a superconductor. The most complicated 
aspect of the phase diagram of the noncompact model is 
the direct phase transition between the NOR phase and 
the FSC phase, which is considerably more difficult to 
characterize than the 2 = 1 order-disorder transition in 
the dissipative 3D XY model. 

A (2+l)-dimensional model with bond dissipation has 
recently been considered as an effective theory of quan- 
tum criticality at optimal doping in high-T c cupratesP^ 
The claim of this work is that the phase correlators of 
the model at the critical point decay algebraically as 1/r 
while they are short-ranged in space. Such a phase tran- 
sition would be an example of local quantum criticality. 
Monte Carlo simulations on the (2+l)-dimcnsional quan- 
tum rotor model gives an order-disorder transition in the 
3D XY universality class, which is quite different from 
local quantum criticality. From the results of the present 
paper, it appears to be important to specify whether the 
phase variables are compact or noncompact, c.f. Fig. |13| 
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The phase transitions separating the CSC phase from 
the FSC phase, or the CSC phase from the NOR phase, 
are not of the type described in Ref. |T7} To verify 
whether or not the remaining phase transition separat- 
ing the FSC phase from the NOR phase is an example of 
local quantum criticality one would ideally need a single 
order parameter measuring spatial and temporal corre- 
lations in phases, 6. Since we do not have this at our 
disposal, we have not been able to determine what sort 
of universality class the critical line separating FSC and 
NOR belongs to, apart form concluding that it is not in 
the 2D XY or 3D XY universality class. However, the 
spatial correlation functions presented in Fig. |T2| suggest 
that the NOR-FSC transition line is not a line with local 
spatial phase correlations. 

We end with an important remark on the temporal 
phase fluctuations we have focused on in this paper. The 
quantity W^ g in Eq. 6 measures temporal fluctuation 
in phase gradients Ad, aefined on a spatial bond of the 
lattice. One could also study a corresponding measure 
of temporal fluctuations of the phases 9 themselves. We 
have done this, and find the following. In all parts of the 
phase diagram in Fig. [T] the quantity 

W2(N r ) = ±(£{6 r --8) a \. (17) 

diverges with N T . This underlines that the instantons 
or instanton-like objects we have discussed in this paper 
are temporal fluctuations in phase gradients A9, not in- 
stantons in phases 9. On the other hand, the hclicity 
modulus Eq. ^measures long-range or quasi-long-range 
spatial ordering of phases 9, and we find such orderings 
in the FSC and CSC phases. Thus, the FSC phase does 
not exhibit 3D XY ordering. It features spatial ordering 
of 9 and A9, but temporal ordering only of AO. This 
supports the statement made above, that the NOR-FSC 
transition is not in the 3D XY universality class. It is a 
new type of phase transition involving a complicated in- 
terplay between spatial point- vortices and instanton-like 
excitations in A9. 



VIII. CONCLUSIONS 



are highly fluctuating in imaginary time and the system 
therefore exhibits metallic behavior. Increasing the dis- 
sipation strength drives the system to a phase transition 
where the phase differences A9 are localized into one of 
the minima of the Josephson potential. This localiza- 
tion of bond variables in imaginary time occurs simul- 
taneously with an onset of rigidity towards phase-twists 
across the spatially extended system. We identify this 
phase with a fully bond-ordered superconducting state 
(FSC). 

At strong coupling and weak dissipation we identify an 
intriguing phase exhibiting finite phase stiffness and al- 
gebraically decaying spatial correlations. The imaginary- 
time direction remains disordered with wildly fluctuating 
bond differences. This dimensionally reduced phase is 
referred to as a critical superconducting (CSC) phase. 
The finite hclicity modulus in this phase indicates that 
the system may sustain a dissipationless current going 
through the entire JJA. There are, however, voltage fluc- 
tuations present which in principle should make it ex- 
perimentally distinguishable from a fully bond-ordered 
superconducting phase, and also distinct from the more 
standard 3D XY ordered fully superconducting state 
where even the phases 9 are ordered in all directions. 

We have found no signs of a phase which is temporally 
ordered (in the sense of having a bounded W^ g ) and pro- 
liferated vortices. Such a phase would naturally facilitate 
the observation of local quantum criticality in which a 
spatially disordered and temporally (quasi-)ordered sys- 
tem disorders in the imaginary-time direction. 
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The model discussed in this paper describes a two- 
dimensional array of quantum dissipative Josephson 
junctions. By extensive Monte Carlo simulations we have 
shown that the model features three distinct phases (see 
Fig. [TJ featuring different behaviors of spatio-temporal 
fluctuations. We have quantified these fluctuations by 
the mean square fluctuation W\ e , Eq. Q6J) , and the spa- 
tial hclicity modulus T r , Eq. (.4). 

The normal phase (NOR) is found at weak dissipation 
and weak Josephson coupling strength. In this phase, 
the spatial hclicity modulus is zero, signaling a vanish- 
ing stiffness to infinitesimal phase twists on each Trotter 
slice. The phase differences of the individual junctions 



Appendix A: Reformulating the noncompact degrees 
of freedom 

To gain further insight into the three phases reported 
in this work and the transitions between them, we con- 
sider the following decomposition of the phase degrees of 
freedom: 

0x,r -> 0x,r + 2nk x . T . (Al) 

The noncompact starting point is thereby exchanged 
for a compact phase field, 9 G [— tt,tt), plus an integer- 
valued field, k, keeping track of the specific 2ir interval 
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the original variable belonged to. In the partition func- 
tion, this reformulation amounts to 



Z=IV9e' s = Y[(d6^ T )e- s (A2) 
E / = E f II ( d ^K*- 

r ; .i J rr. n •'—IT „ „. 



Note that fc is defined on every point in space-time and 
has nothing to do with the winding number found in 
some realizations of quantum rotor models with compact 
phases. 

It should also be noted that the 27r-periodic spatial 
interaction is only sensitive to the 9 field. Also, the com- 
pactness of 9 enables the identification of vortices in this 
field in a similar way as discussed in connection with the 
classical 2D XY-model, Eq. 0. The finite T x observed 
in the CSC and FSC phases may thereby be attributed 
to phase coherence in 9. In addition to the vortex degrees 
of freedom found in the classical version of the system, 
the noncompactness of the quantum version introduces 
an additional degree of freedom (fc) associated with the 
tunneling of bond variables from one minimum of the 
extended Josephson potential to another. 

In the NOR phase, we found T x — 0, which may be 
understood as a phase featuring proliferated vortices of 
the 9 field, as well as proliferated instantons in A6. In- 
creasing the Josephson coupling (for small a) drives the 
system into the CSC phase with T x ^ 0, which corre- 
sponds to a binding of vortices into dipolcs. Nonethe- 
less, the bond variables remain anomalously diffusive, 
W\ e oc lnAV, in both the NOR and the CSC phase. At 
strong coupling, the bond variables tend to stay in the 
vicinity of the minima of the potential wells, A9 ss 27rAfc. 
From the viewpoint of the reformulated variables, the de- 
localized bond variables in the CSC phase is an expres- 
sion of an unbroken symmetry Afc — > Afc + A, where 
A is an integer. Moreover, the integer field fc may be di- 
rectly connected with the instanton charges in the strong 
coupling limit by Afc T+ i — Afc T = nj. The derealiza- 
tion of A9 manifests itself as proliferated instanton/anti- 
instanton configurations in this regime, and the CSC 
phase may therefore be described as an instanton gas. 
This is illustrated in the topmost curve of Fig. [2] illus- 
trating the quantum paths of A9 in the CSC phase of 
Fig. [IJ 

For weak Josephson coupling, the excitations in the 
imaginary-time path of A9 are strictly speaking not well 
described by topological instanton defects. This is quite 
evident from the topmost curve of Fig. [8] describing 
the quantum paths of A9 in the NOR phase of Fig. [T| 
Nevertheless, as the quantum fluctuations still respect 
the same symmetry Afc — > Afc + A, we choose to refer to 
such excitations as instantons also in the NOR phase. In 
the FSC phase, on the other hand, localization of bond 
variables implies that the symmetry is broken for both 
weak and strong coupling. Starting at large K and large 
a, the picture is therefore as follows: In the FSC phase, 



both the defects associated with 9 (vortices) and with 
fc (instantons) are absent or tightly bound. Lowering a 
into the CSC phase, instantons are proliferated while the 
vortices remain bound. Lowering K from the CSC phase 
into the NOR phase, the vortices proliferate as well. 



Appendix B: Correlation functions in the NOR 
phase 



It has recently been propose d 19 * 25 ! that the metallic 
state of Josephson junction arrays might exhibit non- 
trivial behavior. Here, it was argued that the (0+1)- 
dimensional constituents of the array may slide pas t each 
other in what was denoted a "floating phase "P3 Simi- 
lar dimensionally decoupled phases are also believed to 
be relevant to other physical systems such as layered 
superconductors^ and stacks of two-dimensional arrays 
of membrane proteins. 

These papers employed a renormalization group anal- 
ysis to show that the spatial coupling between the su- 
perconducting islands is perturbatively irrelevant on the 
disordered side of the transition. They also calculated 
the correlation functions Eq. ([8]) and (9| in this regime 
and found that they had a form that indicated unconven- 
tional, purely local fluctuations. Monte Carlo studies^of 
a single resistively shunted Josephson junction also indi- 
cated that a similar form of correlation functions could be 
found in (0+l)D systems as well. The correlation func- 
tions employed in these analysis featured a noninteger 
parameter q that was introduced to probe fluctuations 
with another periodicity than the underlying Josephson 
potential. In the presence of a finite Josephson potential, 
expectation values such as (e iAe *<*'- T ) will generally not 
be equal to zero in any phase. This is, however, due to 
the corresponding symmetry being explicitly - and not 
spontaneously - broken, and has consequently nothing to 
do with a phase transition. The parameter q was there- 
fore introduced to assure correlation functions decaying 
to zero in the disordered phase. Similar correlation func- 
tions have also been considered before in investigations of 
roughening transitions of crystal surfaces with quenched 
bulk disorder.^ We will refer to them as fractional cor- 
relation functions. 

Fig. 14 shows both spatial and temporal correlation 
functions, Eq. ([8]), at a dissipation strength deep in the 
NOR phase where the Josephson potential is expected 
to be irrelevant) 2 ^ and we are far away from the phase 
transition at a — Oii\ The correlation functions in the 
bottom row include a noninteger factor q = 1/3, the top 
row shows the correlation functions without (q = 1) this 
noninteger factor. Comparing the correlation function of 
the temporal direction with the spatial direction for q = 
1/3, it is clear that the spatial and temporal behaviors of 
the system appear completely decoupled.^ As we discuss 
below, this local behavior of the fractional correlation 
functions is misleading. 

App. [A] introduced a reformulation of the phase vari- 
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FIG. 14: Correlation functions, Eq. (fSp , in both space (left 
column) and imaginary time (right column) in the normal 
phase. The relevant coupling parameters are a = 0.012, K = 
0.4 and K T = 0.1. There is a pronounced difference between 
the correlations along the spatial and temporal directions for 
q = 1/3. 



ables that clarifies the difficulties concerning the con- 
struction of a globally order parameter describing our 
system. The reformulation of the phase variables also of- 
fers an alternative viewpoint on the fractional correlation 
functions. For example, imagine a 2D XY model, Eq. 
0, being formulated with noncompact phase variables 
instead of the standard compact variables. In the par- 
tition function, the summation over k is trivial, yielding 
only a renormalization of the ground state energy, be- 
cause there is no coupling between different k sectors in 
the action. The remaining integration over 9 is the parti- 
tion function of the ordinary 2D XY model. When per- 
forming Monte Carlo simulations on the 2D XY model 
with a noncompact formulation of the phases, we find 
the usual QLRO phase at strong Josephson coupling, in 
which the correlation function Gg(x; q = 1) of Eq. ^ de- 
cays algebraically. However, consider probing the QLRO 
phase with a fractional correlation function, q < 1. This 
correlation function involves contributions from several 
k sectors, ultimately averaging the correlator to zero for 
all distances. The same result holds for the disordered 
phase, and so, although the QLRO phase is phase co- 
herent and the disordered phase is not, the fractional 
correlation function essentially cannot tell them apart. 

Applying exactly the same arguments as above to our 
CSC phase with spatial QLRO, one realizes that the spa- 
tial fractional correlation function will vanish also here. 
In analogy with the classical 2D XY model, we argue 
that this should not be regarded as a signature of com- 
pletely spatial decoupling in neither the CSC phase nor 
the NOR phase. The apparent locality of the normal 
phase, and by extension the corresponding floating phase 
of Ref 1251 is consequently not a result of the dissipativc 
interaction per se. Rather, a floating phase with such 
vanishing spatial fractional correlations follows as a di- 
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FIG. 15: The distribution of A6 T - A0 O , P(A0,t = 125), 
in arbitrary units for a system with K T = 0.1, K = 0.4 and 
dissipation strength a = 0.011. The distribution is extracted 
from the Monte Carlo simulations and is conjectured to follow 
Eq. pT|. 



rect result of the noncompactness of the phase variables, 
which in turn is caused by their coupling to a dissipative 
bath. 

In the following we provide supplementary details re- 
garding the fractional correlation functions. To be spe- 
cific, we will investigate the fractional bond correlation 
functions Gas(t;<z) more carefully, and prove that a 
power-law tail is expected in the weak dissipation regime. 
We first consider the distribution function P(A9,t), as 
was also the starting point of Ref. EH This function 
describes the diffusion of the phase difference A9 T with 
respect to its value at r = 0. The distribution broadens 
for increasing r and is illustrated for a fixed imaginary- 
time distance in Fig. [15] We find that the distribution 
function can be very well fitted by the functional form 



P{A9, t) = P c 



(Bl) 



where Pq is a normalization constant. The distribution 
is made up of a sequence of sub-gaussians with standard 
deviation a centered around the minima of the Josephson 
potential. In addition, there is an overall gaussian con- 
volution characterized by a standard deviation uq. We 
find empirically that whereas a is dependent on K, it is 
independent of the distance r in imaginary time. The 
overall variance Q{t) of the distribution, as defined by 



G(r) = ((A9 T - A9 ) 2 ) 



(B2) 



grows logarithmically with r. This variance can further- 
more to a very good approximation be identified with the 
variance <7q of the convolution function. 

The calculations in Ref. [HD were based on a strong- 
coupling limit for the distribution function, with an ad- 
ditional assumption that the spatial coupling will renor- 
malize to zero regardless of its bare value. For large val- 
ues of K, we have demonstrated that the system will 
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eventually reach a superconducting state, i.e., the CSC 
phase or the FSC phase, for all a > 0. It is also clear 
from Fig. [15] that there is an appreciable broadening of 
the sub-gaussians (<r > 0) compared to the delta-function 
distribution implicit in the strong-coupling limit (er — > 0). 

We next consider the implications of a finite a on the 
correlation function G/±Q{r]q). Assuming Eq. (Bll, we 
calculate 



j q (Ae T -Ae )\ 



- ^(7 K,q 



E 



(B3) 



where k = a^/(a^+a 2 ). The sum over n has been traded 
for an integral at the cost of introducing an integer Pois- 
son summation variable s. The n- variable is subsequently 
integrated out. Comparing with Eq. (12) in Ref. [HE the 
broadening of the sub-gaussians has introduced an overall 
prefactor and a multiplicative adjustment of the Poisson 
summation variable. The strong-coupling result is eas- 
ily recovered in the limit a — > 0. In the limit t — > oo, 
the term with the slowest decay is dominant, hence the 
sum may be substituted by the term with the smallest 
(q— s/k) 2 . For a logarithmica lly d iverging oq, we also 
have k — y 1, meaning that Eq. (B3l is a scale- free power 
law in this limit. 

In Fig. 16 we show a plot of G?a0(t;q = 3/4) and 
the two terms from Eq. (B3l corresponding to s = 
and s = l. We have set Oq equal to Q(t) as measured 
from the Monte Carlo simulations in order to compare the 
analytical result Eq. (B3l with the fractional correlation 
function found numerically, a is specified from fitting 
Eq. (Bl) to data from Monte Carlo simulations. At 
short distances the s = term is still contributing, but a 
clear cross-over to the dominant s = 1 term is visible for 
larger values of r. The excellent fit between the curves 



validates the functional form of the distribution (Bl ) 



It is interesting to compare the behavior presented 
above with available numerical results for a single re- 
sistively shunted Josephson junction. Ref. |TT| reports 
temporal fractional correlation functions in a (0+l)D sys- 
tem that are power law in much the same way as those 
in Ref. UHl They also report a logarithmically diverg- 
ing MSD, but only at the phase boundary. This is in 
contrast to the results presented in Sees. IVI] and |IV| 



where we find logarithmic growth as a generic feature 
of the weak-dissipation phases. Following Ref. 11, it 
is natural to consider the possibility that a logarithmi- 
cally diverging MSD is the signature of critical behavior 
for models describing Josephson junctions. A logarith- 
mically diverging MSD follows from a logarithmically di- 
verging Q(j), and we have shown that the latter generates 
fractional bond correlators that are algebraically decay- 
ing in imaginary time. A possible scenario could be that 
the increased dimensionality of the problem has damped 
the fluctuations such that, in contrast to the single junc- 
tion, the entire weak-dissipation regime features critical 
temporal correlations of the bond variables. However, 
we expect such a critical phase to produce divergent sus- 



10 



< 



1(T 



10 




-s = 
-s = 1 

-G Ae (T;q = 3/A) 



FIG. 16: (Color online) The unequal-time bond correlation 
function, Gab(t; q = 3/4), for a system with K = 0.4, K T = 
0.1, a = 0.011, N = 20 and iV T = 250. The black curve is 
the correlation function Eq. Q sampled directly from the 
Monte Carlo simulations. The red (lowermost gray) and blue 
(upp ermost gray) curve are the s = and s = 1 terms of Eq. 
(|B3[), respectively, and are calculated as explained in the text. 



ceptibilities of the action. The simulations do not sup- 
port this scenario and we find nonanalytic xs only at the 
points a = ai 1 ^'^. Thus, a power-law form of the tem- 
poral fractional bond correlators can not necessarily be 
ascribed to critical behavior of the system. 
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